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Abstract. We review the basic concepts of all-order calculations in Quantum 
Chromodynamics (QCD) and their application to collider phenomenology. We start 
by discussing the factorization properties of QCD amplitudes and cross-sections 
in the soft and collinear limits and their resulting all-order exponentiation. We 
then discuss several applications of this formalism to observables which are of great 
interest at particle colliders. In this context, we describe the all-order resummation 
of event-shape distributions, as well as observables that probe the internal structure 
of hadronic jets. 


1. Introduction 

The CERN Large Hadron Collider (LHC) has recently resumed its operations after 
the first long shutdown. Run I of the LHC was extremely successful, with the 
milestone discovery of the long-sought-after Higgs boson by the ATLAS [1] and CMS [2] 
collaborations. Moreover, many other aspects of the Standard Model have been probed 
in a previously unexplored energy regime. As a result, after LHC Run I, the Standard 
Model appears now as a fully consistent and highly-successful theory of particle physics. 
However, very strong evidence of its incomplete nature already exists, above all the 
fact that no Standard Model particle appears to be a good candidate for dark matter. 

This situation results into a two-fold task for LHC Run H. On the one hand, 
precise measurements of the Higgs boson properties are necessary, in order to verify 
whether it is fully responsible for electroweak symmetry breaking as predicted by the 
Brout-Englert-Higgs mechanism [3-6]. On the other hand, searches for new particles 
and possible inconsistencies in Standard Model predictions are going to be pushed to 
a new energy frontier. In order for these tasks to be successful and to fully exploit 
the LHC physics potential, new and more accurate theoretical and experimental tools 
have been, and are being, developed. The aim of this topical review is to discuss 
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one of these theoretical tools, namely all-order calculations in perturbative Quantum 
Chromodynamics (QCD), the theory of strong interactions. 

The LHC collides protons, which are strongly interacting, and further strongly- 
interacting particles are abundantly produced in every such collision. Therefore, 
careful studies of QCD radiation in Higgs production and new physics processes 
can be exploited in order to better understand their properties. Moreover, the 
possibility of making discoveries at the LHC depends on our ability to separate new 
and rare phenomena from an overwhelming background, which is often several orders 
of magnitude bigger than the signal. This background consists of Standard Model 
processes and its dominant component comes from strong interactions. It follows 
that precision QCD is mandatory at the LHC. Furthermore, we need an accurate 
understanding of physics in the presence of disparate energy scales, which range from 
the unprecedentedly large colliding energy, through the electroweak scale, all the way 
down to hadron masses. The appearance of multiple scales renders perturbative QCD 
calculations unreliable at any hnite order. Therefore, an all-order re-organization of 
the perturbative expansion is necessary in order obtain reliable theoretical predictions. 

Finally, despite being often taken for granted, it is far from trivial that precise 
perturbative QCD predictions can be compared to observables which are measured in 
terms of hnal-state hadrons. As we are going to discuss in the following, there exist 
special classes of observables for which corrections due to the hadronization process 
scale with inverse powers of the energy involved in the hard scattering and can therefore 
be considered reasonably under control at sufficiently high energies. Moreover, since 
it is not possible to compute these corrections from hrst principles, the use accurate 
parton-level predictions when comparing to data, allows us to constrain different ways 
of modeling the hadronization process. We will briefly come back to this in section 4. 

In this topical review we discuss the basic concepts behind all-order calculations 
in QCD. More than one technique to perform such computations has been developed, 
and here we focus on analyses which directly make use of QCD matrix elements in the 
relevant soft and collinear limits. Analogous results can also be obtained using the 
methods of Soft-Collinear Effective Theory (SCET) [7-13] (see Ref. [14] for a recent, 
extensive, review). Moreover, recent studies have discussed the equivalence of the two 
methods [15-20]. 

Even if we limit ourselves to resummation in direct QCD, where factorization 
is established for amplitudes and cross-sections in the soft or collinear limit, we still 
encounter different ways of obtaining all-order results. Generally speaking, we think 
one can identify an American school and a European (including Russia) one, although 
clear distinctions are difficult and sometimes misleading. The former, typically 
describes resummation by introducing non-local correlation operators, such as Wilson 
lines, and exploits their renormalization group evolution. The latter instead resorts 
to a more iterative procedure, directly identifying factorization and exponentiation 
properties of QCD matrix elements and cross-sections. In our presentation we will 
mostly follow this second approach. 

In the hrst part of this review, we are going to discuss basic properties 
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of gauge-theory amplitudes in the soft and collinear limits, specihcally Quantum 
Electrodynamics (QED) and QCD, starting at one loop in section 2 and then moving 
to an all-order analysis in section 3. In the second part of this review, we are going to 
focus on phenomenological applications. We will start with an analysis of event-shape 
variables in section 4, discussing thrust in some detail. We will then generalize our 
discussion to non-global observables in section 5. We will hnish with a selection of 
more advanced topics in jet physics and jet substructure in section 6, and a discussion 
about the limitations of factorization theorems upon which resummation is founded 
in section 7, before concluding. 

For each topic, we are going to describe in some detail results that are well 
established in the literature, providing an extensive list of references. We are also going 
to provide short summaries of recent developments. We hope this way to stimulate 
the curiosity of both beginner and expert Readers. 

2. Infrared divergences 

In order to better understand the origin of the corrections that we wish to resum, we 
start our discussion by making some considerations about the properties of scattering 
amplitudes in the soft and collinear limits. In this regime, matrix elements involving 
massless particles may exhibit divergences. More specifically, the presence of soft 
divergences is related to the emission or exchange of particles with vanishing four- 
momentum, and it is associated with the presence of massless vector bosons. Such 
divergences occur even when the matter particles are massive. Collinear divergences 
are instead related to the splitting of particles at small angles, and are strictly present 
only when all involved particles are massless. Although these statements are rather 
general, in the following we will focus on QCD with some digression to QED, where 
several aspects can be simplihed because of its Abelian nature. 

An important point to make is that the strong coupling grows in the infrared 
and therefore perturbative QCD loses its predictive power at long distances. However, 
even without worrying about non-perturbative contributions, we also have the issue of 
how to define asymptotic states in gauge theories (QCD and QED alike), because of 
the presence of soft bremsstrahlung. 

Our strategy is to define measurable quantities that are not affected by infrared 
and collinear (IRC) divergences. The Bloch-Nordsieck (BN) [21] and Kinoshita-Lee- 
Nauenberg (KLN) [22,23] theorems state that observable transition probabilities are 
free of IRC singularities. In our analysis we are going to consider safe observables, i.e. 
measurable quantities that do not spoil the above theorems. We will come back to a 
more precise definition of IRC safety in sections 2.2 and 4. It is worth pointing out 
that from an experimental viewpoint, the finite resolution of the detectors acts as a 
regulator, thus preventing the occurrence of actual singularities. However, this in turn 
would be reflected on a possibly strong dependence of theoretical predictions on the 
detector resolution parameters, which we wish to avoid. 

Even if we focus on IRC safe observables, from a practical viewpoint, when 
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Figure 1. Feynman diagrams contributing to the cross-section of e“'"e —>■ qqg. 


computing Feynman diagrams, we need a recipe in order to deal with potentially 
divergent intermediate contributions. We adopt the following 

(i) sum over all possible and indistinguishable processes that lead to the same hnal 
state conhguration; 

(ii) compute cross-sections with an infrared regulator (e.g. in dimensional 
regularization with d = 4 — 2e dimensions, where e < 0); 

(iii) select observables for which the 4-dimensional limit is hnite, i.e. that do not spoil 
BN and KLN theorems; 

(iv) interpret these results as perturbative estimate of the corresponding hadronic 
cross-section measured in experiments. 

Thus, if we limit ourselves to the aforementioned class of safe observables, 
perturbation theory will provide us with a hnite result. However, as we shall discuss 
below, the cancellation of IRC singularities can leave behind a hnite, but potentially 
large contribution, making the perturbative expansion unreliable at any hnite order. 
Our task is to identify these problematic IRC contributions and to re-organize the 
perturbative expansion in such a way that they are accounted for to all orders. Even 
though the structure of IRC singularities in QCD is universal, in this section we will 
consider a simple example to illustrate their structure, namely the cross-section for 
electron-positron (e'''e“) annihilation into hadrons. 
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2.1. Hadronic cross-section in e'^e annihilation 

The general expression for the cross-section e’''e“ —>■ hadrons can be written as 

^had ('S) <^had,0 T ^2^) ^had,l T ^2^) *^had,2 T ^ (®s) i (^) 

2 

where ttg = is the strong coupling constant. The leading order (LO) contribution 
(Figure 1.1) can be cast as 

” “ ( 2 ) 


dTTCr \ 

O'had,0 = ' 


3s 


q ) 


where we have introduced the hne structure constant a, the quark electric charge Cg, 
the number of colors Nc and the squared center-of-mass energy s. At next-to-leading 
order (NLO) the real and virtual corrections have similar Feynman diagrams with 
different kinematics (Figure 1.3-1. 4 and 1.2 respectively). The real contribution is 
proportional to: 

■S 13 S 23 / 


E 

Spin 


IM 3 


oc 


■S 13 ^ _l_ 

S 23 Si3 


where Sij = (fcj -|- and Sijk = {ki -|- kj -|- k^)'^. The final state quark and anti-quark 
have momenta ki and /c 2 , whereas the gluon has momentum k^. The sum in front of 
the squared 3-particle hnal state amplitude M 3 indicates the sum over the hnal state 
particle spins and the average over the initial once. The amplitude is singular in the 
limits: 

S 13 —t 0 , S 23 —t 0 -V^ i ?3 —)■ 0 , 013 —)■ 0 , 023 —)■ 0. (4) 

Rewriting it in terms of energy fractions Xi = 2ki ■ Q/Q^^ where Q = fci -|- A :2 + ^3 is 
the center-of-mass energy, we can explicitly split the amplitude into a singular and a 
non-singular part 


Spin 


iMsl oc 
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2. (5) 
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The above equation makes all singular limits explicit: the hrst term encodes the soft 
singularity and correspond to the Altarelli-Parisi splitting kernel 

1 -1- (1 - xf 


Pgq 
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which describes the probability of a collinear splitting. We can therefore write 

,2 1 . . 1 


E 
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IM 3 


OC 


1 — Xi 


Pgq ( 3 ^ 3 ) + 


1 - 0:2 


Pgq {xfi) - 2. 


( 6 ) 


( 7 ) 


The hrst two terms describe the two possible collinear limits. The last term is instead 
hnite. In order to evaluate the phase-space integral, we use dimensional regularization 
in d = 4 — 2e dimensions. The singularities now appear in terms of poles in e:| 

Cp (^ 


<ad,l = t^had.O 


F(1-£) 


2 3 2 19 . 

- + - - TT^ + — + O £ 


( 8 ) 


f A step-by-step derivation of this result can be found e.g. in Ref. [24]. 
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The virtual amplitude in Feynman gauge receives a contribution only from the vertex 
correction, which, once squared and integrated over phase space, gives 


virtual 

^had,l 


Cp f 47rp^\ ^ / 2 

r(i^V W V ^ 



8 + O (£) 


(9) 


Therefore, the NLO contribution to the total cross-section is 


_ _ _real i —virtual _ 

<^had,l — <^had,l + ‘^had,l ~ <^had,0 


c. 


r(i-T) 




2 + o(>;) 


( 10 ) 


and the total cross-section itself can be written as 

O'had = O'had.O ^ ("^s)) ) (H) 

which is hnite. The real and the virtual contributions have the same singularity 
structure with opposite sign and the IRC poles in e cancel leaving a finite result, 
as stated by the KLN theorem. This cancellation occurs because in the soft and 
collinear limits both the real and the virtual amplitudes become proportional to the 
Born amplitude. In what follows we are going to investigate this last point in more 
detail. 


2.2. Factorization in the soft limit 


In this section we study the factorization properties of real and virtual amplitudes in 
the soft limit. We still focus on the above example and we start considering just one 
of the two diagrams (the one depicted in Fig. 1.4) contributing to the emission of a 
gluon from the quark anti-quark pair (for this reason the matrix element is primed). 
In order to stress the special role of the gluon in the soft limit, we are going to indicate 
its (soft) momentum by q, rather than k^. We have 


M's = tl (fci) (g) 


+ 


g—>0 


(g kif' + is 

h 


t1g,g,^u{ki)Y£,.{q) . , , ,2 , • 

[q + ki) 

K 


M 2 

M 2 


= t‘lg,ia^-^—ef,{q) M 2 , 


( 12 ) 


where the factor kf/{ki ■ q) is called eikonal factor and is the color matrix. We 
have also used fairly standard notation for the Dirac spinor u{k) and for the gluon 
polarization vector s^{q). In the last step the Dirac spinor was absorbed in the 2- 
parton matrix element M 2 and therefore we dropped the tilde on it. For the full 
amplitude we find 


M3 {pi,V2 ki,k2,q) 



(g) Sf, (g) M 2 {pi,P2 -t ki, /cs), 


(13) 


where we have introduced the eikonal current 

ktl 
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n=l 


■ q 


( 14 ) 
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It is important to note that the factorization does not depend on the internal structure 
of the amplitude. From the physical point of view, this reflects the fact that the 
large wavelength of the soft radiation cannot resolve the details of the short distance 
interactions. Squaring the amplitude, we obtain the following factorization property 
for the soft real emission: 




E ta ,, 


h . h 




{km ■ q){kn ■ q) 
{ki ■ fcs) 


{ki ■ q){k 2 ■ q) ' 

The soft approximation can be applied also to the loop amplitude. In this limit we 
can in general neglect powers of the loop momentum q in the numerator if 
furthermore in the denominator we can use the fact that <C ki-q. The loop correction 
to quark-antiquark pair production is therefore proportional to 

d'^g u (ki) Y {^1 + <0 Y [i - h) (^ 2 ) 


I = 




{27iY [(g + ki)‘^ + ie] [(g - ^ 2 )^ + is] [g^ + is] 
Yq {ki ■ ^ 2 ) [u {ki) Y'v (^ 2 )] 


(16) 


(27r)‘^ [q ■ ki + ie] [-q ■ k 2 + is] [g^ ie] ’ 

The result in d = 4 dimensions can be evaluated in the center-of-mass system of ki 
and k 2 , with 

= Ei(l,0,0,1) , = ^2 (1,0, 0,-1) , g'^ = (go,g) with g = (gl,g^) , (17) 

where q± is the vectorial transverse loop momentum and we dehne gj_ = | gl |. 

We thus obtain 

d^g 2dgo[h(/ci)7^u(fc2)] 


I = gi^Ypi-z) 


(2vr)4 (go -q^ + ie) (-go - g^ + ie) (g^ 


Qz 


gi + ie) 


(18) 


Eq. (18) has four poles in the complex go plane at 


qo = qz-is, qo = -q^ + ie, qo = ± {\q\ + ie) . (19) 

Closing the contour from below we hnd 

^ [h (fci) 7^^^ (^ 2 )] j 

where the second integral is a pure phase 

J ( 27 r)^ {qz - ^s) (gi) J ql + e^ J {fipf g± J ( 27 r)^ g± 

We are going to refer to the above contribution as Coulomb phase. We note that 
the above phase always cancels when considering physical cross-section in Abelian 
theories like QED. However, it can have a measurable effect in QCD cross-sections, in 
the presence of a high enough number of harder colored legs, as discussed in section 4.2. 


- {ki ■ k2) 


2 |g I (/ci ■ g) (/c 2 ■ g) i.qz - is) {qf) _ 


,( 20 ) 
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Collecting real and virtnal contribntions together, we can compnte the NLO 
distribntion of an observable v by introdncing an appropriate observable-fnnction 

u^{{k}y. 

^ “ oTi / ^$2 1^2!^ (^ 1 ; ^2) + 

2 r d^q 1 2 ^ ih ■ k2 

J (2vr)3 2|g| ^ {h ■ q) {h ■ q) 

X [U3{ki,k2,q)-U2{ki,k2)] . (22) 

From the last eqnality we can derive the following important conclnsions: 

• for a complete cancellation of the IRC contribntions it is important that the 
observable is infrared and collinear safe, i.e. according to the dehnition of Ref. [25], 
that it satishes: 

Um+i (• • • 5 k{^ kj ,...) y Um (• • • 5 d" kj, ■ • •) if ki II kj^ (^3) 

Ura+l (• • • ) • • •) ^ Um (..., ki-i, ...) if fcj —0. (24) 

These limits have to hold not only for a single particle, bnt for an ensamble of 

partons becoming soft and/or collinear. IRC safe properties of jet cross-sections 

and related variables, snch as event shapes and energy correlation fnnctions were 
hrst stndied in Refs. [26-28]. We note here that there exists a wealth of observables 
that are of great interest despite them being IRC nnsafe. Generally speaking, these 
observables reqnire the introdnction of non-pertnrbative fnnctions to describe 
their soft and/or collinear behavior. For example, lepton-hadron and hadron- 
hadron cross-sections are written as a momentnm-fraction convolntion of partonic 
cross-sections and parton distribntion fnnctions. Arbitrary collinear emissions 
change the valne of the momentnm fraction that enters the hard scattering, 
resnlting in nn-cancelled collinear singnlarities. Finite cross-sections are then 
obtained by a renormalization procednre of the parton densities. Similar sitnations 
are also enconntered in hnal-state evolntion, if one is interested in measnring 
a particular type of hadron (see e.g. [29]) or if the measurement only involves 
charged particles [30,31]. Furthermore, we mention that recent work [32-34] has 
introduced the concept of Sudakov safety, which enables to extend the reach of 
(resummed) perturbation theory beyond the IRC domain. 

• in the case of inclusive observables, for which Um {ki,. ■ ■, km) = 1 for all m, the 
cancellation is complete. Consequently, the total cross-section remains unchanged 
by the emission of soft particles, as it should. 

• in case of an exclusive (but IRC safe) measurement, although the singularities 
cancel, the kinematic dependence of the observable can cause an unbalance 
between real and virtual contributions, which manifests itself with the appearance 
of potentially large logarithmic corrections to any orders in perturbation theory. 

These large logarithmic contributions spoil the perturbative expansion in the strong 
coupling and must be resummed to all orders in order to obtain reliable theoretical 
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predictions for exclusive measurements. In the following, we will discuss how the all¬ 
order behavior of QCD in the IRC regime can be systematically captured, enabling 
us to obtain perturbative predictions in kinematical regions where the fixed-order 
expansion in the strong coupling breaks down. 


3. A first look at all-order resummation: the coherent branching 
algorithm in the leading collinear approximation 

In this section we present a first all-order analysis of gauge theories (QED and QCD) 
based on the formalism of generating functionals [35]. Our discussion follows the 
presentation given in the lecture series of Ref. [36]. 


3.1. Soft emission probability in QED 

We start by considering the factorization of a matrix element in the presence of a soft 
emission as derived in Eqs. (13) and (14), but this time we restrict ourselves to the 
case of an Abelian gauge theory like QED. Labeling the hard momenta with {fcj} and 
the momentum which becomes unresolved with q, and considering the possibility of 
emitting a photon by all the charged external legs, we can write 

M^+i {{ki} , q) gjf^ (g) (g) ({A:*}) , (25) 

where is the eikonal current 

ur- 

( 26 ) 

i 

The sign of the lepton charge e* is assigned depending on whether the momentum ki 
is incoming or outgoing. Because of charge conservation we have 

= (27) 

i 

Squaring Eq. (25) and taking into account also the phase space we obtain an expression 
for the factorization of the cross-section 


dffm+i {{ki} , q) = {{ki})y. dWi (q) , 


(28) 


where dWi is the single photon emission probability. It can be written as 

dlVl(,) = [dg]J»‘(,)^^„„J■'(,)9^ 


(29) 


where is the sum over the photon polarizations. Because of charge conservation 
only the term proportional to the metric tensor is different from zero when d^iy is 
contracted with the eikonal currents. The result of this contraction can be cast into 
(we choose here an explicit reference frame for the momenta and define u to be the 
energy of the soft photon) 


(g) Ju (?) 


2 1 — cos Oij 

^2 (l-cosB,,) ' 


(30) 
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Note that Eq. (30) contains also the effects of the interference between emissions from 
different external lines. The physical conseqnences of this interference is important, 
we therefore analyze the last eqnation a bit further§. We consider the term in the snm 

_ (31) 

(1 — COS dig) (1 — COS 0jq) ’ 


and split it into two terms 




m 


where 


1E“ = - 

2 


+ 1-T 

1 — cos 




i^j 


(32) 


(33) 


^ik f COS 

The two terms on the right hand side of Eq. (32) have angnlar ordering properties. 
If we write the angnlar integration in terms of the polar and azimnthal angle with 
respect to the momentnm ki and we integrate over the azimnth, we find 


-271 




1 


= 
27r 


1 —COS 6iq 
0 


if diq < Oj^j , 


(34) 

otherwise. 

This means that each of the two terms describe radiation which is confined in a cone. 
The physical interpretation behind this is that photons at larger angles are nnable to 
resolve the pair of partons (?, j) as separate charges. We consider now again the single 
photon emission probability of Eq. (29) in the leading collinear approximation of q 
being collinear to kj. This means that we consider the angle 9jq to be much smaller 
than any other relative angle. In this limit we can recast Eq. (30) as 

(1 — cos 9ij) 


j2 ! \ ^ 79-^0 1 

J (?) = 


2ej 


(1 — cos 9jq) 


E 


e* 


(1 - cos9iq) 


+ 0 ( 1 ) 


(35) 


The terms of O (1) are regular in the limit 9jq —?• 0 and can be neglected whereas, in 
the considered limit, the sum can be approximated as 
(1 —cos 6'ij) 






j- 


(1 - cos9ik) 

Writing the differential phase space element as 

dfiq a;^da;d(;/)dcos6' 


[dq] = 


{27rf 2u (27r)^ 2u 


we can express the single photon emission probability as 


dW'i (9) = - E 

TT f ^ 


, den d9l 


u 




p© (^'max — 9jq) , 


(36) 


(37) 


(38) 


where we have taken the limit of small angle 9, as appropriate for the collinear region; 
^'max is of the order of the other angles 9ij. This means that we can approximate the 
single photon emission probability by considering a sum over radiation emitted by 
independent sources, where the destructive interference, which cancels the radiation 
at large angles, is approximated by the angular ordering constraint. 

§ The following considerations are derived in detail also in Chap 5.5 of the book by Ellis, Stirling 
and Webber [37]. 
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We can now generalize the emission of a single soft photon to multiple emissions from 
the same fermionic line. Since photons do not carry electric charge and since the 
emissions we consider are soft, they leave both the charge and the momentum of the 
emitting particle unchanged. Therefore we can write the multiple emission probability 
as 

1 "" 

dW„ ~ - TT dWi (g,) , (39) 

n! 

i=l 

where the prefactor is the symmetry factor for n identical bosons. 

The emission of multiple soft photons to all orders can conveniently be described 
by introducing a generating functional: 

OO « 

[u (g)] = 1 + X] / dlTn (gi,..., gn) M (gi) -...■« (gn), (40) 

n=l 


where for each photon we introduced an arbitrary weight, which acts as phase-space 
constraint similarly to the observable function Um introduced in the previous section, 
but in the soft-collinear limit. From Eq. (40) one can recover any emission probability 
by successive differentiation at m = 0: 


dW (gi,... ,gn) 


5$ 


6u (gi) ■ ... ■ 6u (gn) 


1 


Using Eq. (39) we can rewrite Eq. (40) as 


(41) 


OO -| n r p 
n=l i=l L"' 


dWi (qi) u (gi) 
= exp <; / dfUi (g) u (g) 


n=l L-" 


— / dWi (qi) u (qi) 


(42) 


This equation gives the corrections due to real soft photon emission to a squared 
matrix element in QED. Form the last line we observe that it is simply given by 
exponentiating the lowest order contribution. This was derived for the first time in 
the 1960s in Ref. [38]. To hnd the total correction due to soft emission we also have 
to consider the virtual contribution. In the previous section we explicitly computed 
the loop contribution and found that the total (unconstrained) soft emission has a 
vanishing effect. We can therefore exploit this result by imposing: 


<1>[“(4)1L., = 1- («) 

This is often referred to as the unitarity condition and allows us to correctly normalize 
<h, which can hnally be written as 

^real 


$[M(g)] = 


,|,real ^ 


= exp 


dWi (g) [^(g) - 1] 0 (Q - w) 0 {u9 - Qq) 


(44) 
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Figure 2. Schematic representation of soft emissions in QED and QCD. On the 
left independent emission from different external legs in QED. On the right angular 
ordered emissions along an external leg in QED (top) and in QCD (bottom). 


In the last equation we have introduced two constraints in terms of step-functions, 
the hrst one gives an upper limit on the energy considered in the soft approximation 
Q, which is typically of the order of the hard scale considered in the process. The 
second 0-function introduces an arbitrary lower cutoff Qo for the photon transverse 
momentum. The dependence on Qo does drop off when IRC safe observables are 
considered. Note that by requiring that no real radiation is emitted, i.e. setting 
M (g) = 0, we obtain the so-called the Sudakov form factor 

$ [u = 0] = A (Q, Qo) = exp J dhhi (g) 0 (Q - w) 0 {coO - Qo) |, (45) 

which describes the probability that no emission takes place and is fully determined by 
the virtual contributions. In the limit of vanishing lower cutoff Qo, the exponentiated 
integral diverges and the Sudakov form factor vanishes as well. This means that the 
probability of not observing any radiation is zero. In other words, in any scattering 
process there must be some radiation, which can however be arbitrarily soft. 

When considering multiple emissions from different charged fermions we have to 
remember the effect of interference, which can be approximated by the constraint on 
the angular ordering. For this case the generating functional can be written as 

[u (g)] = JJ [Q, 0max; u (g)] , (46) 

i 

where for each charged fermion external line we have 

[Q,0max;M(g)] = exp j-y —J —[u{q)-l]Q{ue-Qo)^. (47) 
A schematic representation of the angular ordered emissions is given in Fig. 2. 
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We now extend the generating functional method to QCD. The non-Abelian structure 
of the color charge leads to some complications in comparison with the derivation of the 
previous section. Charges have a matrix representation in color space and therefore the 
factorization in the soft limit is more involved. The example case used in section 2.2 
to derive the factorization in the soft limit has a very simple color structure. For 
processes with more colored particles Eqs. (14) and (15) have to be generalized. The 
eikonal current can be written as 

. ki-q 

where the Tj’s are abstract color operators (see for instance Refs. [39-43]), which 
satisfy color conservation 

^T,|M^({fc})) =0. (49) 


In the last equation | Mm{{k})) is a vector in color space. The square of an operator 
Tj gives the Casimir of the SU{Nc) representation 


Cp if i is a quark or anti-quark, 
Ca if i is a gluon. 


(60) 


while, once a color basis is hxed, the color products Tj ■ are represented by 
matrices, which are, in general, non-diagonal. We note that in the large-A'c limit, 
often considered, for instance, in Monte Carlo parton showers, the off-diagonal entries 
of these matrices vanish. 

The factorization of an (m -|- l)-parton matrix element in the soft limit can thus 
be written as 

’^■>3 

An important consequence of the non-Abelian structure is that factorization is 
generally incomplete because of the presence of color matrices. Moreover, we also 
have to consider the fact that soft radiation can come from hard gluons and that, no 
matter how soft the emitted gluons are, they will always carry away charge. The exact 
pattern of soft gluon radiation is therefore more involved than in QED, and in order 
to keep this presentation simple, we limit ourselves to the leading collinear behavior, 
where we have remarkable simplihcations. We shall come back to the more general 
case in section 4, where we consider the resummation of a specihc class of observables. 

In the leading collinear approximation, the probability for a single emission is very 
similar to the expression we have previously derived for QED, Eq. (38), 

dITl (g) = - V Q —^0 (0max - e,Q , (52) 

vr ^ u 

I ‘y 


where Ci = Cp yn. case parton Ms a quark or an anti-quark and C* = Ca in case parton 
i is a gluon. As in the case of photon emission interference effects can be captured 
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by means of angular ordering. This is strictly true only at the logarithmic order we 
are considering. As we shall see later, wide-angle soft radiation, whose contribution is 
subleading, makes this pattern much more involved. 

As already anticipated, when considering multiple gluon emission color 
correlations cannot be neglected. After the hrst splitting the total color charge will 
be shared among the two partons and further radiation can be emitted from either 
of them. This more complicated radiation pattern can be simplihed with the help of 
color coherence. In fact, as derived in Eq. (34) above, soft radiation cannot resolve 
the details of the interaction which happens at shorter distance and higher momentum 
scale. Therefore a soft gluon emitted at an angle 6i will only “see” the total color 
charge of the radiation emitted at smaller angles 6 < 6i, which corresponds to the 
charge of the hardest parton involved [44-46] . Iterating this argument we can describe 
multiple soft QCD radiation in the leading collinear approximation in a similar way 
to the QED case, remembering that all the radiation needs to be angular ordered. 
Furthermore, each emitted parton can act as a further emitter, as represented on the 
lower right picture of Fig. 2. The hnal equation for the generating functional in leading 
collinear approximation in QCD can thus be written as 

['^^(?)] ~ ^ki [A'fcij ^max) (^)] j (53) 


where however is it not possible to write a closed expression for <hfc., but only an 
iterative one [47-49] 


[Ek,,0i^^;u{q)] = exp 




s'^kj 

TT 


dujq 

j max 



UJq 

L 

01, 


^kiq] u] 

-1)} 


(54) 


The hrst term in the second line describes further real radiation which is softer and 
occurs at a smaller angle, the second term encodes the virtual corrections and can be 
derived via a unitarity argument as done in the previous section. Eqs. (53) and (54) 
constitute the essence of the so called coherent branching algorithm [50,51]. They are 
the starting point for constructing Monte Carlo parton shower algorithms [52] such as 
Herwig [53] and for the all-order resummation of the leading logarithms, which can 
be systematically improved, as we will discuss in the following. 

From Eq. (54), we can easily rederive the single emission probability, which reads 


dPFi = 




4i=0 


5u{k) 


= <l>fc [Ffc, 0] 


71 


dCU d0? 

) Jo ^kq 


TV 


dUn 


UJn 


d9l^ 


Olq 


^ ^kq-) 0) 


A(St,9^,,0) (55) 


In the last equation the hrst and the last Sudakov form factors describe the probability 
of evolving without branching of the mother particle before and after the emission. 
The emission probability is given by the term in square brackets, which includes the 
probability for the emitted particle to evolve without further branchings. The double 
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emission probability can be obtained computing a further derivative. The recursiveness 
of Eq. (54) will generate two terms in that case: one describing an angular ordered 
double emission from the same mother particle and another one describing a further 
splitting of the first emitted gluon. 

The logarithmic accuracy of this description can be straightforwardly improved 
to include hard collinear radiation by considering the Altarelli-Parisi functions, rather 
than their soft limit, as splitting kernels, while the inclusion of soft emissions away 
from the collinear limit is more involved. Both of these issues will be discussed in 
the next section, where we are going to present all-order results for a wide class of 
observables, i.e. event shapes. 

4. Event shapes in e+e“ collisions 

Thus far we have discussed general properties of QCD matrix elements in the soft and 
collinear regime, showing how perturbative matrix elements factorize in those limits, 
leading to the possibility of all-order calculations in QCD. We now turn our attention 
to the all-order behavior of specihc classes of observables. 

We would like to start with an example that although simple, contains most 
of the ingredients which form the building blocks of resummation technology. A 
good candidate for this program is actually a fairly large class of observables called 
event shapes, which aim to measure the energy flow of an event. In order to avoid 
possible complications due to initial-state radiation and the role of parton distribution 
functions, we begin by considering event shapes in annihilation. Not only this 

simplifies our analysis, but it also has phenomenological value. Although event shapes 
can be defined and studied at hadron colliders, e.g. [54,55], they are best measured at 
lepton colliders, such as LEP, where they are used to determine the strong coupling 
constant, e.g. [56-62]. The absence of some non-perturbative phenomena, such as the 
underlying event and pile-up, makes event shapes at lepton colliders an invaluable tool 
to study the energy-momentum flow of strongly interacting final states. 

We are going to consider observables V = V {pi,... ,Pn), which are functions of 
the hnal-state momenta pi,... ,Pn, and that satisfy two basics properties: 

(i) infrared and collinear safety: an event shape is IRC safe if its value v does 
not change in the presence of an arbitrary number of strictly soft or collinear 
emissions (see section 2). 

(ii) globalness: an observable is defined to be global if it is sensitive to radiation 
anywhere in phase-space. 

As already discussed in section 2, the first condition ensures cancellation of soft and 
collinear singularities between real and virtual corrections, so that the observable v 
can be computed in (hxed-order) perturbation theory. Moreover, it is also possible to 
show that non-perturbative corrections to event-shape distributions, such as the ones 
due to the hadronization process, are suppressed by inverse powers of the hard scale 
(see e.g. Refs. [63-70]). Many observables of phenomenological interest, especially 
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at hadron colliders do not obey the second condition and they are therefore termed 
non-global [71]. Their all-order structure is more intricate and will be discussed in 
section 5. 

To ensure the feasibility of resummation, the two conditions above have to be 
made more restrictive, as specihed in Refs. [72-76]. In fact, one needs to impose: 

(i) recursive infrared and collinear safety: an event shape is rlRC safe if two 
conditions are met 

(a) in the presence of multiple soft and/or collinear emissions, the observable 
should have the same scaling properties as with just one of them; 

(b) for sufficiently small values V = v oi the observable, there exists some £ -C 1 
such that e can be chosen independently of v and emissions below ev do not 
signihcantly contribute to the observable. 

(ii) continuons globalness: an observable is defined to be continuously global if 
its scaling with respect to the transverse momentum of a soft and/or collinear 
emission is the same everywhere in phase-space. 

We refer the interested Reader to Ref. [75] for a detailed discussion. 

In our analysis, it will prove very convenient to consider the cumulative 
distribution for an event shape, defined as integral of the differential distribution up 
to the value V = v, normalized by the total cross-section: 

(56) 

<^o Jo 

In the next section we consider an explicit example for such an event shape, namely 
thrust T [77], defined in Eq. (58). Many more event shape observables exist and we 
refer to [75,76,78,79] for an exhaustive list. 

We adopt the convention that the logarithmic accuracy of a calculation is 
determined with respect to the logarithm of S, which after resummation can be written 
in exponentiated form as 

S(n) oc exp {Lgi (usL) -F g 2 {a^L) + Usgs {asL) ...} , (57) 

where L = In^. A leading logarithmic (LL) calculation then determines g\, which 
resums contributions in In E to all orders in perturbation theory. In this 

review we will be mostly considering next-to-leading logarithmic (NLL) accuracy, 
which corresponds to computing gi and g 2 , i.e. resumming all a^ln*^ ^ contributions 
to liiE. Thus the large parameter which is resummed is a^L < 1. For completeness 
we remind that also other conventions can be adopted. In particular, for observables 
which do not exponentiate, it is possible to reorganize the perturbative expansion 
in terms of powers of a^L'^. The hrst term in this convention, also known as tower 
expansion [80], sums all double-logarithmic terms the second one the terms 

asL{asL‘^Y~^, and so on, assuming that the large parameter is OsT^ < 1. 
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Figure 3. Typical thrust values for back-to-back two-jet configurations (left) and 
three-jet hnal states (right). 


4 . 1 . The thrust distribution 


In order to make this exposition as pedagogical as possible, while still discnssing 
an observable of phenomenological relevance, we begin by considering as a concrete 
example the event shape thrust [77]. The NLL resummation of thrust, and other 
related event shapes, was first performed in Ref. [51], which we closely follow in this 
presentation. 

Given a collection of hnal-state momenta fipi] = we dehne 


T = max 


Ei \Pi 


■ n\ 


Ei \Pi\ 


= max 


Ei \Pi 


■ n\ 


Q 


(58) 


where the second equality holds if all particles are massless, i.e. fif = H and Q is 
the center-of-mass energy of the colliding electrons. The three-vector n = ht that 
maximizes the sum in Eq. (58) is called the thrust axis. As depicted in Fig. 3, thrust 
measures how uniform radiation is distributed in the event, with T ~ 1 indicating an 
event which is two-jet like. For convenience, the variable r = 1 —T is often introduced. 

Before discussing the resummation of the thrust distribution, let us analyze its 
kinematics in the presence of soft and/or collinear radiation. Let us call the plane 
orthogonal to the thrust axis This plane divides the event into two hemispheres 
Si and 82 - We also dehne hemisphere momenta: 


qi= '^Pi = zip + qti + zip 

ieSi 

q2= X] Pi = Z2P + qt2 + Z2p, (59) 

ieS 2 

where p,p are lightlike momenta. We now state two properties of the thrust axis 
Ht [51], namely that no pi lies in and that the hemisphere three-momenta 
are aligned with ut, i.e. qn = qt 2 = 0. Note that g/ 7 ^ 0. Thanks to these 
properties, and using energy-momentum conservation together with the above Sudakov 
parametrization we have 


r = 1 - — (|gi ■ utI -b \q2 ■ nr|) = 1 - 
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= + A + 'M'] = , ‘^P-Eies,Pi <ih^\ 

\Q^) Q 2 + Q 2 ^ ) 

n 

= J2^'{k,). ( 60 ) 

i=l 

where in the last line the sum runs over n soft and/or collinear emissions, which we 
have labelled with momenta ki. Note that we have neglected recoil effects, which 
are beyond NLL accuracy [51,74,81]. Thus, we have managed to turn the task of 
resumming the thrust distribution, into the resummation of the invariant masses of 
the two hemispheres. Moreover, the second and third last lines of Eq. (60) show 
that to the accuracy we are working at, the hemisphere mass and thrust are additive 
observables, i.e. the contribution of n emission is the sum of the contributions of each 
emissions. 

We can now turn our attention to the actual resummation of the thrust 
distribution. We start by considering the real-emission contribution to the cumulative 
distribution Eq. (56), in the appropriate collinear limit. At NLL accuracy it can 
be written as a product of independent angular ordered emissions, as given by the 
generating functional derived in the previous section: 




oo 


n=0 


n\ 


n 

i=l ■ 




kl 


2 n 


F^^Pgq[Zi)Q 




i=l 


( 61 ) 


where the transverse momenta ku are measured with respect to the thrust axis. We 
have introduced the splitting function Pgq{z) = and the n\ accounts for ordering 

of the emitted gluons. The strong coupling Og is evaluated in the so called CMW 
scheme [50] 

«s = af® 1 K = Ca “ y) “ 


i.e. it accounts inclusively for secondary branchings at NLL level. Note that K is the 
coefficient of the two-loop cusp anomalous dimension||. Factorization properties of 
Eq. (61) appear to be spoiled by the presence of the ©-function that constrains the 
measured value of thrust and prevents us from summing the series into an exponential 
factor. As it is often the case in this kind of situation, the kinematic constraint can 
be diagonalized by considering a conjugate space. In particular, we introduce the 
following integral representation: 


0 



du 

2Tiiv 


n 



2 = 1 


( 63 ) 


|j The dynamics of soft radiation can be computed considering correlators of Wilson lines with a cusp. 
Their renormalization leads to the cusp anomalous dimension, which controls the renormalization 
group evolution of these amplitudes [82-84]. In the CMW scheme the first subleading contribution 
associated with soft emissions is accounted for with a redefinition of the strong coupling, allowing to 
capture the full NLL dependence. 












QCD resummation for hadronic final states 


19 


where the dependence oi V on p and p is understood. Thus, by inserting the above 
expression into Eq. (61), we are able to perform the sum over the number of emissions 
and we obtain 


= 


dz/ 

2nii> 

dz/ 

2nii> 




n\ 

n =0 2=1 




kl 27r 


27r 


l/T 

e exp 




2C, 


afikt 


kf 271 ^ 271 


-Pggi^ 


-uV(k) 


(64) 


We now have to consider virtual corrections. In the soft or collinear limit, these 
are equal and opposite to the real-emission one. However, no observable constraint 
is present because the kinematics is fixed to Born level, leading straightforwardly 
to exponentiation. Thus, putting together real and virtual corrections, the all-order 
cumulative distribution reads 

dz/ 


E(r) = 


27rzz/ 


VT 

e exp 


d^ 


dk^ d(j) 

k? 271 


2C,^fi>p„{z) - 


1 ) 


(65) 


We now want to compute the resummed exponent to a fixed logarithmic accuracy. In 
this example, and for most of this review, we are working at NLL accuracy. To this 
purpose, we evaluate the integrals over the gluon phase-space, making the following 
approximat ions: 

(i) we evaluate the running coupling with the two-loop QCD fi function; 

(ii) we also note that the factor — l) essentially acts as a restriction on the 

allowed phase space. To NLL accuracy we can replace it with [85] 

_ 1 ) ^ _0 {y{k) - e-^^z/-^) , 


where is the Euler constant. The generalization of this procedure to higher 
logarithmic accuracy is discussed in Appendix A of Ref. [ 86 ]; 

(iii) we integrate the finite part of the splitting function down to z = 0, i.e. 


B, = I 


We obtain 


E(r) = 


-e exp 

27rzz/ 


— 9l{^v) + 92{^v) 




( 66 ) 


with Xy = CKg/^olnz/ and Og = Q;g(Q). The function gi resums LL contributions, while 
g 2 NLL ones: 


^i(A) = - 

92{X) = - 


Cf 

CpK 


[(1 - 2Xy) ln(l - 2Xy) - 2(1 - A,) ln(l - A,)], 


[21n(l-A,) -ln(l-2A,)] + 


CpB, 

TTfio 


ln(l - Xy) 


Cpfii 

TTfil 


- -hT(l -2A 
2 ^ 


21n(l — Xy) — ln(l — 2Ai.) -|- ln^(l — Xy) 
CplE 


- 2 - 


tt/^o 


[ln(l - Xy) - ln(l - 2Ai.)]. 


(67) 
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In order to perform the inverse Mellin transform, we expand the exponent in powers 
of In z/ around In z/ = — In r 



( 68 ) 


with A = asfioL and L = — In r and the dots stay for higher derivatives, which give rise 
to sub-leading contributions. Thus, we have reached the main result of this section: 


Eq. (68) represents the NLL resummed cumulative distribution for thrust. The result 


can be upgraded to sometimes called NLL' accuracy by considering the 0{as) constant 
contributions, schematically S —)■ (1 -|- asgo^i)T,. 

Some comments about this result are in order. Eq. (68) is expected to capture the 


dominant physical effects in the region where the logarithms are large, i.e. ~ 1, 


but the use of perturbative QCD is still justihed. Note that the functions gi have 


logarithmic branch cuts starting from A = 1/2, which corresponds to r = e ^“=^ 0 . This 
singularity originates from the Landau pole of the QCD running coupling and heralds 
the breakdown of the perturbative approach. Even before reaching this singularity, as 
the value of r decreases, it becomes sensitive to emissions with transverse momentum 
in the non-perturbative region, and non-perturbative power corrections must be taken 
into account. A thorough discussion of these effects goes beyond the purpose of 
this review and we refer the interested reader to the original literature for more 
details [64,87-91]. Let us just briefly mention that several approaches exist. Their 
basic common idea is to refrain from doing the substitution (ii) described above in 
Eq. (65), and instead to split the integral in the exponential into a perturbative and a 
non-perturbative part. In the so-called dispersive model [63-65] the non-perturbative 
part is treated by dehning an effective coupling, which is supposed to be hnite in the 
infrared region and which can be written in terms of a dispersive relation. Recently 
this model was extended to match NNLL-I-NNLO accuracy for thrust [62]. Other 
alternatives are the single dressed gluon approximation [90,92,93], which assumes the 
existence of a reordering of the perturbative series in a so-called skeleton expansion, 
or the introduction of a shape function [89, 94, 95], which is a non-perturbative 
function which admits an operator dehnition in terms Wilson lines, and was recently 
extended and used in the context of SCET [60,91]. The inclusion of non-perturbative 
effects is important in the context of strong coupling determinations, where also their 
universality and the effect of hnite mass corrections can be tested [66,69,96-98]. 

The resummed result on its own is also inadequate for describing the region where 
r is an order one quantity, because the soft and/or collinear approximation on which 
resummed calculations are based upon breaks down. This is the domain of validity of 
hxed-order calculations, which can be matched to the resummation in order to obtain 
solid theoretical predictions over a vast range of r. As an example, in Fig. 4, we show 
the matched distribution for thrust. The calculation is actually performed at a higher 
accuracy than the one presented above: the resummation is performed to NNLL and 
it is matched to a NNLO hxed-order calculation [81]. A similar result with part of the 
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Figure 4. The distribution of the thrust event shape at NNLL matched to NNLO, 
in red, compared to the fixed-order (NNLO) result, in green. The bands indicate the 
theoretical uncertainty. Plot taken from Ref. [81]. 


NNNLL contribution has been computed also in SCET [59]. 

It is also useful to discuss the limitations of the framework in which the above 
calculations has been carried out, as well as possible extensions. Firstly, the original 
derivation of the all-order thrust distribution was based on the coherent branching 
algorithm described in section 3, which exploits color-coherence of azimuthally- 
integrated matrix elements that are relevant for one- and two-jet observables [50,51]. 
A framework that enables one to perform resummed calculations in the presence 
of many hard legs is desirable but not straightforward beyond LL accuracy. While 
single-logarithmic terms arising from hard collinear emissions are still captured by 
a probabilistic method, the correct description of large-angle soft radiation requires 
to go beyond it. This issue was addressed by different groups. For instance, NLL 
resummation of near-to-planar 3-jet events was performed in Ref. [99], while general 
frameworks to resum soft large-angle logarithms in processes with many legs was 
worked out in Refs. [100-103] and Ref [104], The hnal result is actually remarkably 
simple and it amounts to consider soft emissions from each of the dipole that constitute 
the system of n hard partons, weighted by an operator describing the color exchange. 
We will come back to this point in the next section. 

The second limitation that we wish to stress has to do with the nature of the 
observable that we wish to resum. Central to the above analysis was the identihcation 
of the integral transform to diagonalize the observable Eq. (63), which relied upon 
the additive nature of thrust. Different, and possibly more complicated observables, 
require multiple integral transforms, restricting the set of resummable observables to 
the ones for which a theorist’s ingenuity has not failed to provide the correct integral 
transform. This is not a pleasant situation and a more automated method is highly 
desirable. Such a framework was devised in Refs [72-76] and will be the subject of the 
next section. 
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Caesar [72-76] is a framework (and a computer program) that allows one to perform 
the resummation of a large class of observables, namely global event shapes, to NLL 
accuracy. We begin this discussion by considering, as in the previous section, processes 
which at Born level feature two hard massless partons (plus eventually color singlets, 
e.g. photons, Higgs or electroweak bosons) and we denote the set of Born momenta 
with {p} We consider positive-dehnite observables V that measure the difference 
in the energy-momentum flow of an event with respect to the Born conhguration. In 
particular, we focus on observables that vanish when evaluated on Born conhgurations, 
V ({p}) = 0. The previously discussed variable r = 1 — T, which vanishes in the 2-jet 
limit, is an example of such observables. 

We consider the cumulative distribution for the event shape v in the presence of 
n soft or collinear emissions ki 

S(n) = 

°° 1 r 

e (« - r({fc})) . (69) 

n=0 i 

Henceforth, the dependence of V upon the Born momenta is understood, i.e. V ({fc^}) = 
V{{p},{ki}) The hrst term corresponds to virtual corrections to the Born process, 
which exponentiate in the soft or collinear limit, as previously discussed. is 

the matrix element squared for the real emission of n soft or collinear partons off 
the two hard partons. The key idea is to further divide the real-emission phase- 
space by introducing a resolution parameter £. Resolved emissions give V > ev, 
while unresolved ones V < ev. Because of rlRC safety, unresolved emissions do not 
signihcantly contribute to the measured value of the event shape and therefore, like 
virtual corrections, they exponentiate 

S(n) = [ Yl[dki\M^{ki), (70) 

n=0 i 

where we have introduced the shorthand notation 

r \{[dki]M\ki) = [ l[[dh]M^{k,)e {v - V{{t^)) 

Jev ^ J ^ 

xQ{V{{h})-ev) (71) 

It is now natural to combine virtual and unresolved contributions: 

S(n) = (72) 

where we have introduced 

jr(y) = (73) 

n=0 ” i=l 

The hrst contribution only depends on one emission and, once evaluated to a hxed- 
logarithmic accuracy, it is the generalization of the resummed exponent in the 


^ In principle one should consider the set of momenta {p} after recoil, but this effect is beyond NLL. 
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resummation of thrust Eq. (68). It has the following physical interpretation: given 
an emission that sets the value v of the event shape, it vetoes further emissions which 
would contribute to the event shape more than v. The second term instead takes into 
account the effect of multiple emissions that equally contribute to the event shape. 

We now briefly describe how to evaluate both the resummed exponent and the 
function iF to NLL. Consider a single emission with momentum k, which is soft and 
collinear to a hard leg 1. We can parametrize the value of the event shape in the 
presence of this emission as follows 

v(k) = d,(^^'^ pl'>) , (74) 


where k^ \ and denote transverse momentum, rapidity and azimuth of the 
emission, all measured with respect to parton 1 . Q is the hard scale of the process 
which can be set, for instance, equal to the partonic centre of mass energy, i.e. Q = i/s. 
Note that the requirement of continuous globalness requires all the a coefficients to be 
the same, while IRC safety imposes a > 0. 

Thus far we have considered, mainly for simplicity, only the case of an 
underlying Born with two hard partons. However, as we are about to evaluate the 
resummed exponent, we can lift this restriction and consider an ensemble of m hard 
partons. The main complication we have to deal with is the presence of many color 
conhgurations, which renders the virtual correction matrices in color space. Therefore 
the exponentiation in Eq. (72) has to be understood in a formal way. The calculation 
of F instead does not change at this logarithmic accuracy. We have ^ 

E(u) = -(Mo|e-^'e-^|Mo)^, (75) 

<^0 


where \Mq) is the Born amplitude, i.e. cxo = {Mo|Mo). The real part of the resummed 
exponent has a structure similar to the one we have encountered for thrust 


E E (-T. 

dipoles ij legs l£ij 


2 fc? 2ir 2ir ' ' ' ' 


X e 







(76) 


1 (0 

where rj = ^\n and p is the appropriate reduced splitting function. The integral 
above corresponds to the shaded area in Fig. 5. Note that we have arbitrarily separated 
the integration ai p = 0, identifying a region where emissions are collinear to leg i and 
leg j respectively. The precise position of this boundary is beyond NLL accuracy. 
Moreover, we have also made use of the color product Tj ■ Tj operators introduced in 
section 3.3 and applied in the context of resummation for the first time in Ref. [104]. 

The resummed exponent R has also an imaginary part known as Coulomb phase 
(see discussion after Eq. (18)) which is proportional to Tj-T^, where i,j are both final- 
state partons (or both initial-state partons in the case of hadron collision). This phase 


Beyond NLL accuracy considered here, one needs to introduce path-order exponentials in order to 
deal with the exponentiation of non-commuting matrices. 



QCD resummation for hadronic final states 


24 



Figure 5. The emission phase-space in the ( 77 , In At/Q) plane, as parametrized by 
the Caesar formula Eq. (74.) We note that, in the soft and collinear limit, emissions 
are uniformly distributed in this plane. Figure taken from Ref. [74]. 


cancels in the prodnct of exponentials when m < 4, while gives a physical contribntion 
to the distribntions in processes with fonr or more hard legs. 

It is now relatively straightforward to compnte the one-emission integrals; nsing 
color-conservation Eq. (49) we are able to cast the Caesar resnmmed formnla Eq. (72) 
in the same form as Eq. (68). The LL fnnction gi that appears in the resnmmed 
exponent generalizes to 


gfiasL) = ^gifiasL) = - ^ 


C, 


1=1 


— {a + bi — 2A) In I 1 — 


27r/3^6; L 
2A 


[a — 2A) In ( 1-'j 


(77) 


a bi ^ 

The above LL resnlt consists of a snm over all hard partons, the dependence on the 
color is trivial and only enters throngh the Casimir of each leg I, {Cp for a qnark leg, 
Ca for a glnon leg). The thrnst resnlt is recovered by considering only one qq dipole 
and by setting a = bi = digi = 1, I = 1,2. The resnlt for the NLL fnnction g 2 has a 
richer strnctnre: 
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The hrst term in the sqnare brackets in Eq. (78) contains the two-loop contribntion 
to the DGLAP splitting fnnction in the soft limit and to the QCD /^-fnnction 
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a bi 



(79) 


The second term in the square brackets in Eq. (78) instead captures hard collinear 
emissions to a quark leg {Bg = — |) or to a gluon leg {Bg = —nfio); we have also 
introduced 


T(L) 

Di 


nfio 1 - 2asfioL' 



(80) 

(81) 


where Ei is the energy of leg 1. The contribution S instead captures the effect of 
large-angle soft emission. At NLL, it is the only contribution for which we have to 
keep track of the matrix structure in color space described above. In order to explicitly 
evaluate this contribution, one needs to fix a color basis, compute the representation 
of the color products Tj ■ and the color-decomposed Born amplitude. Although 
straightforward, these steps can be tedious and laborious because the dimensionality 
of the basis quickly increases with the number of hard legs. For this reason, until very 
recently, only calculations with up to four hard legs were performed *. 

We now turn our attention to J-, which describes the effect of multiple emissions. 
This contribution only starts at single-logarithmic level, therefore, in order to perform 
a NLL calculation, we can consider all the emissions in Eq. (73) to be soft and 
collinear and widely separated in rapidity. Moreover, the parametrization of the 
observable simplifies as well, because we can ignore the overall normalization digi{(j)). 
These considerations lead to a simpler expression in the case of event-shape variables 
considered here (see Ref. [74] and Ref. [109] for details): 


n=0 i=l I 

v{{p},m) 


X 0 1 — lim 


D-5-0 V 


(82) 


where Q = Vi/v is the fractional contribution to the observable due to emission i. The 
above expression contains only NLL terms and it is suitable for numerical evaluation 
via a Monte Carlo generator. 

Thanks to the Caesar framework, it has been possible to perform resummed 
calculations for a variety of event shapes, for which no explicit integral transform to 
diagonalize the observable was known. Moreover, as we shall see in the next section, 
this framework can be extended, with minor modihcations, to event shapes in hadron 
collisions. 

Finally, we note that in some cases further singularities can appear inside the 
physical region of the phase space. These singularities originate for example when 
the phase-space boundary for a given number of partons lies inside the one for a 
larger number of partons. In this case, if the observable is discontinuous in that point. 


* Refs. [105-108] considered the resummation of 2 —3 scattering processes, in the limit where one 
of the final-state partons was a soft gluon. 
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Figure 6. Distribution of the directly global thrust minor in pp collisions computed 
in three different approximations. The fixed-order results are computed at LO (in 
green), and NLO (in blue). The resummed and matched distribution is instead 
computed with Caesar at NLL-I-NLO (in red). Plot taken from Ref. [76]. 


logarithms will appear which have to be resummed to all orders to make the observable 
hnite there. One speaks in this case of a Sudakov shoulder. These logarithms can 
also be resummed [110], and in the very general case of cuts applied in experimental 
analyses, this task is usually fulhlled by parton-shower Monte Carlo generators. 


4-3. Extension to lepton-hadron and hadron-hadron collisions 


The above analysis for event shapes in lepton colliders can be extended to the case 
of lepton-hadron and hadron-hadron collisions (even though the framework we are 
discussing may break down at sufficiently high orders, see section 7). When computing 
cross-sections with one or two hadrons in the initial state, parton distribution functions 
(PDFs) are introduced in order to describe the long-distance non-perturbative 
physics of the initial state. PDFs obey the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equation [111-113], which resums single logarithms due to the 
emissions of collinear partons off the incoming legs up to a factorization scale, which 
is usually taken to be of the order of the hard scale of the process Q. 

However, the phase-space restriction that is imposed by measuring an event shape 
prevents us from inclusively integrating collinear emissions above kt ~ Qv^+'^i , because 
such emissions would result in too high a value for the event shape. Therefore, we have 
to correct the PDF scale choice in the Born cross-section (which only knows about Q) 
in order to account for this effect. This leads to a new contribution to the NLL function 
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which effectively replaces the PDF q^^fix, Q^) evaluated at the hard scale, with the one 
evaluated at a lower scale, set by the event shape. We further note that the virtual 
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part of collinear corrections are already accounted for by the Bi term in Eq. (78). The 
actual proof of this result is not entirely straightforward and we refer the interested 
reader to Appendix E of Ref. [74]. 

Phenomenological studies of event shapes at lepton-hadron colliders using 
resummed perturbation theory have been reviewed in Ref. [114] (see Ref. [115] 
for recent developments), while the hadron-hadron case was studied in Ref. [76]. 
An example of resummed and matched result is shown in Fig. 6, where the all¬ 
order (NLL-I-NLO) prediction for an event shape is compared to its LO and NLO 
approximations. The observable chosen as an example is the directly global thrust 
minor as measured in pp collisions. 

A very interesting topic, both theoretically and phenomenologically, which would 
deserve a review on its own, is the resummation of the transverse momentum {Qt) 
distribution of electro-weak hnal states, such as the Higgs boson or a lepton pair 
produced via the Drell-Yan mechanism. 

The literature on Qt resummation is vast and since the seminal paper Ref. [116], 
there has been a continuous effort in producing accurate theoretical predictions that 
can describe the experimental data. For example, high logarithmic accuracy [117-123] 
has been achieved and computer programs that allow one to obtain NNLL-I-NLO 
predictions for the Qt distribution in case of colorless hnal states in hadron collision 
exist, e.g. [118, 119, 124-128]. Very recently hrst Qt resummation results were 
computed also for colored hnal states, and in particular for heavy quarks [129-131]. 

Novel observables such as the ot [132] and [133] variables that exploit angular 
correlations to probe similar physics as Qt, while being measured with much better 
experimental uncertainty, have been introduced. This triggered theoretical studies to 
extend the formalism of Qt resummation to these observables [126,134-137] (for a brief 
review see Ref. [138]). The experimental resolution of cff is so good [139-143] that the 
theoretical uncertainty of the state-of-the-art NNLL-I-NLO calculation is much larger 
than the experimental one, calling for improved theoretical predictions. 

An important issue in the context of Higgs physics is the ability of separating 
events according to their jet multiplicity. In particular, events where the Higgs 
is produced in association with n-jets are usually identihed by vetoing additional 
radiation above a given threshold. Jet vetoes are usually applied to transverse 
momentum variables and therefore the all-order treatment of the 0-jet bin of Higgs 
cross-section is naturally related to the Higgs Qt spectrum itself [144]. This is actually 
true at low-logarithmic accuracy, while other effects, such as the ones due to the parton 
recombination to form a jet, become relevant at higher orders [145-147]. 

Thus far we have discussed the importance of resummation for exclusive processes. 
We also want to mention the fact that there are situations where resummation can 
become relevant even for inclusive cross-sections. Let us consider the ration between 
the typical hard scale of a process over the collider energy, x = Q'^/s. If Q^ ~ s, 
and consequently a: —)• 1, the hnal state is produced near threshold and logarithms 
of 1 — a: become large and may require resummation. As for the aforementioned 
Qt resummation, there exists an extensive literature on threshold resummation, see 
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e.g. [85,148,149]. Here, we simply report the state of the art, which is N^LL for 2 —>■ 1 
processes, such as deep inelastic scattering, Higgs production and Drell-Yan [150-156]. 
NNLL threshold resummation is also included in state-of-the-art theoretical predictions 
for the tt cross-section [157,158]. 

In the opposite, high energy, limit of x —)■ 0, small-x logarithms become dominant. 
Such contributions, which are governed by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
equation [159-164], contaminates both the evolution of the parton densities and of 
the partonic coefficient functions. The resummation of small-x contributions to PDF 
evolution was investigated in the 1990s by more than one group to NLL, see, for 
example. Refs. [165-168] and Refs. [169-174]. The resummation of partonic coefficient 
functions is based on the so-called fcf-factorization theorem [175-182] and it is known 
to LL for an increasing number of cross-sections and distributions [183-188] . Monte 
Carlo programs that incorporate small-x contributions also exist, e.g. Cascade [189] 
and Hej [190]. 

We conclude this brief discussion by mentioning a case where the inclusion of 
all-order effects has proved itself crucial in the context of interpreting a discrepancy 
between theory and experimental data. Both the ATLAS and CMS collaborations have 
reported an excess in the inclusive WW cross-section with respect to the Standard 
Model predictions [191-194]. This has clearly sparked a lot of interests as a possible 
manifestation of new physics. However, careful analyses have shown that the excess 
can be largely explained within the Standard Model, if all-order resummation is taken 
into account when computing the measured fiducial cross-section [195-198]. 

4-4- Recent developments 

Before finishing this chapter, we would like to mention two recent developments 
in the resummation of event shapes. As we discussed in detail, the Caesar 
approach provides a highly-automated framework to compute global event shapes at 
NLL. Moreover, the structure of soft singularities in multi-loop amplitudes has been 
extensively studied in the literature, e.g. [199-202] and their all-order behavior turns 
out to be highly-constraint by symmetries [203-205]. 

However, as we mentioned earlier, in order to correctly capture the logarithmic 
terms originating from emissions of soft gluon at large angles beyond LL, color 
bases have to be identified in order to explicitly compute color operators and Born 
amplitudes. Because of the increasing complexity of these bases, until recent, 
calculations were performed by-hand for relatively low multiplicities, e.g. 2 —)■ 2 and 
2 —)■ 3 QCD scattering. However, color-flow information is built in modern matrix 
element generators, such as CoMix [206] or Madgraph [207]. Recent work has 
solved this problem by constructing and implementing a framework that exploits this 
observation and allows for a calculation of large-angle soft logarithms in an highly 
automated way [208]. The algorithm constructs an appropriate color basis for the 
partonic process at hand, and evaluates color operators and the decomposition of Born 
amplitudes in that basis; color-ordered partial amplitudes are then evaluated using 
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Comix . By merging this recent development together with the Caesar approach, 
NLL resummation in the presence of many hard jets, e.g. 2 —)■ 5 QCD scattering, is 
now automated. 

The second important recent development is the extension of the resummation to 
NNLL accuracy, which was achieved for global 2-jet event shapes in Ref. [109]. The 
resummed distribution can be written as 

E(u) = [J^(A) + (5J^(A)]. (84) 

Thus, in order to achieve NNLL, corrections to both the resolved contribution {6J^) 
and the unresolved one, together with virtual corrections (g^), must be computed. In 
Ref. [109] the calculation of 6iF is performed. We remind the Reader that the NLL 
calculation of iF was essentially done in the soft and collinear limit. This hypothesis 
must be relaxed if we are to achieve NNLL. However, the crucial observation is that 
we are only interested in corrections which are one power of Og higher than NLL. 
Therefore, it is sufficient to relax the hypotheses that went into the NLL calculation 
for one emission at the time. In particular, the ingredients that go into the NNLL 
determination of the resolved emission contributions are [109] 

• exact rapidity bound and running coupling corrections to the soft and collinear 
function F{v) 

• one of the emissions ki is collinear but not soft, generating hard-collinear and 
recoil corrections; 

• one of the emissions is soft but at wide angle; 

• gluon decay is treated non-inclusively, giving rise to a correlated-emission 
correction; 

• two emissions are close in angle (rapidity); this contribution usually vanishes for 
event shapes, but does play a role in other observables, e.g. jet vetoes. 

Correspondingly, NNLL corrections to the resummed exponent, encoded in the 
function must be computed. Although the Authors of Ref. [109] do not provide 
a generic computation for all these contributions, they crucially identify classes of 
event shapes that share the same g^ function, thus enabling them to perform NNLL 
resummation for a fairly large class of event shapes, having as input the relatively few 
observables for which NNLL resummation was already known. 

5. Non-global logarithms 

In the previous section we have presented all-order calculations which allow for 
resummation of large logarithmic corrections to NLL accuracy for a large class of 
event shapes. The two key properties that an observable must satisfy in order for 
that analysis to hold are IRC safety and globalness. In this section, we will relax 
the second hypothesis and we will discuss the resummation of so called non-global 
observables [71,209]. In the Caesar framework, non-globalness can be viewed as a 
particular failure of rlRC safety. Not only this is an interesting theoretical question 
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Figure 7. Kinematic configurations that give rise to non-global logarithms to lowest 
order in perturbation theory. Gluon ki G Si does not contribute to the observable, 
while ^2 G S 2 does. This figure is adapted from Ref. [71]. 


from the point of view of understanding the all-order structure of perturbative QCD, 
but it is also central for phenomenological studies. Many observables studied at lepton 
and hadron colliders are sensitive only to a restricted region of phase-space. For 
instance, when considering the properties of hadronic jets, we are asking questions 
about the energy-momentum flow within a certain phase-space region, the jet, while 
ignoring the outside. We start our discussion with a hxed-order example, which will 
illustrate how a single logarithmic contribution arises in non-global observables, while 
being absent for global ones. 

In the previous section we have analyzed in some detail the kinematics and the 
resummation of thrust. In particular, we have related the event shape r to the sum of 
the invariant masses of the two hemispheres Si and 82 - From that analysis, we have 
learnt that, in order to capture NLL accuracy, it is sufficient to examine the emission of 
soft and collinear partons (real and virtual) as if they were emitted independently. All 
contributions originating from correlated emissions were accounted for by the running 
of the strong coupling. 

Let us analyze the situation in more detail, focussing on the hrst non-trivial order, 
i.e. O (Og). At this order we consider the emission of two gluons with momentum ki 
and k 2 - If the gluons are both emitted off the quark-antiquark line their contribution is 
proportional to the color factor Cfi. We can also encounter a situation where one gluon 
is emitted and then splits into a gluon pair, both ending up in the same hemisphere. 
This contribution is a correlated emission, proportional to CaCf, which is captured 
by the running of the strong coupling in the CMW scheme [50]. Let us also consider 
the case in which the emitted gluon splits into two gluons, each of which ends up in 
different hemispheres. To this order we then have r = -^ + -^+0 ( j . We note that 
the limit in which k 2 is much softer than fci, ^2 gives a vanishing contribution to thrust. 
Therefore, when real and virtual corrections are added together, the cancellation of 
the singularities is complete and no large logarithm appear from this conhguration. 
Thus, because the phase-space for thrust is fully inclusive we are able to conclude that 
correlated emissions are fully accounted for by the running of Og, in the CMW scheme. 

What would happen if instead of thrust, i.e. the sum of the hemisphere masses, 
we were interested in the mass of one hemisphere, say S 2 "! The analysis concerning 
independent emissions and running coupling will go through in a similar way. Indeed, 
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in the previous section, we built the resummation of thrust as the sum of the two 
jet mass contributions. However, the story dramatically changes for conhgurations in 
which gluon 1 is emitted in Si and gluon 2 in S' 2 , as depicted in Fig. 7. In this case 
ki does not contribute to the observable because it is in the wrong hemisphere Si, 
while /c 2 does, regardless of its softness. Cancellation of real and virtual diagrams is 
therefore incomplete and a large logarithm is left behind. 

Let us work through this example explicitly. We need to consider the matrix 
element square for the emission of two soft gluon off a qq dipole, in the limit where 
/c 2 is much softer than ki [210,211]. This can be written as the sum of two pieces: 
independent and correlated emissions 

W = (85) 


where 


ly(ind) 


pp(corr) 


2p■p 2p■p 

p ■ kip ■ kip ■ k2p ■ k2 

2 p ■ p / p ■ ki ^ p ■ ki 
p- kip - ki \p ■ k2ki- k2 p- k2ki- k2 


P-P \ 
p ■ k2p- k2j 


( 86 ) 


(87) 


In order to compute the non global contribution to the cumulative distribution E, we 
have to integrate the above matrix element over the appropriate phase-space. Adding 
together real and virtual contribution, we obtain 

= - CfCa [ [dki] [ [dk2] ^ 2 ) 

Jk^(^S2 Jk2&S2 

xQ[2p-k2- Q^v) = -CfCa y 1. (88) 

Thus, the hemisphere mass distribution receives a non-global contribution starting 
from 0{a‘^), which is alLfi, i.e. it contributes to NLL. We note that the logarithm 
in Eq. (88) originates from the energy integrals: non-global logarithms are related 
to large-angle emissions. The coefficient of this logarithm is hxed by the angular 
integrations, which exhibits an integrable singularity when the gluons’ momenta are 
close together, i.e. at the boundary between the two hemispheres [209]. 

The result in Eq. (88) represents only the leading term at the first order at which 
non-global logarithms appear. In order to achieve NLL accuracy these contributions 
must be resummed to all-orders. As we shall discuss in the next section this is very 
non-trivial, even if our aim is to resum the leading tower of non-global logarithms. 


5.1. Resummation of non-global logarithms in the large Nc limit 

In order to perform an all-order analysis of non-global logarithms, we must consider 
conhgurations of many soft gluons. If we restrict ourselves to considering their leading 
contributions, which we recall is single-logarithmic, we can assume energy-ordering; 
however, no collinear approximation can be made. Thus, we have to describe the 
emission of a softer gluon off an ensemble of large-angle soft gluons. As previously 
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discussed, color correlations make the color algebra highly non-trivial as every emission 
increases the dimensionality of the relevant color space. Moreover, describing the 
geometry of such ensemble becomes also difficult. The approach that was taken in 
the hrst analysis of non-global logarithms was to consider the large-Nc limit. Color 


correlations becomes trivial in this limit because the off-diagonal entries of the color 


matrices vanish. Thus, we are able to write the matrix element square for the n gluon 
ensemble in a factorized way [212] 



p ■ p 


(89) 


p ■ ki^ ki^ -ki^-.-ki^-p 


where the sum is over the n\ permutations. For instance, it is easy to verify that in 
the large Nc limit, Eq. (89) with n = 2 and Eq. (85) coincide. Thus, Eq. (89) leads 
to a simplihed physical picture because an emission off an ensemble of n — 1 gluons 
(plus the two hard patrons) reduces to the emission off each of the n dipoles. When 
the dipole radiates a gluon, it splits into two dipoles, originating conhgurations which 
are determined by the history of the gluon branching. This suggests to make use of a 
Monte Carlo implementation, which enables one to deal numerically with the second 
difficulty we have mentioned, namely the complicated geometry of the multi-gluon 
hnal states. This solution was hrst implemented in Ref. [71] and subsequently used in 
a number of phenomenological applications, e.g. [213-217]. We will come back to the 
numerical impact of non-global logarithms when we discuss the jet mass distribution 
in section 6.1. 

Ref. [218] instead developed a more formal treatment of non-global logarithms. 
Starting from Eq. (89) the Authors of Ref. [218] were able to derive an evolution 
equation, henceforth the Banh-Marchesini-Syme (BMS) equation, which, equivalently 
to the Monte Carlo approach, resums the leading non-global logarithm, in the large- 
Nc limit. In order to write down the BMS equation let us consider an observable n, 
which receives contributions only from emissions within region X, and not from its 
complement X, respectively S 2 and Si in the hemisphere mass example above. For a 
generic pair of primary partons a and b, not necessarily aligned with the thrust axis, 
the distribution of L = — In n obeys the following equation: 



(90) 


We note that the hrst contribution is linear in G and is only sensitive to the emission 
of a soft gluon off the primary dipole. This term generates a resummed exponent 
completely analogous to the one for global events shapes that we have discussed in the 
previous section. However, the striking feature of the BMS equation is its second, non¬ 
linear, contribution, which highlights the non-global nature of the evolution. Despite 
the fact that the BMS equation has been derived a while ago, no closed-form analytic 
solution is known and solutions have been determined either numerically or iteratively. 
For instance, the Authors of Ref. [219] have been recently computed an iterative 
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solution up to five loops by exploiting underlying symmetries of the equation. In 
Ref. [220] the expansion is calculated even one order further by means of a Monte 
Carlo approach. These results have also been conhrmed, and extended to hnite Nc, 
by means of brute force calculations of Feynman diagrams in the soft limit [221]. 

It has been noted [222] that the BMS equation has the same form as the Baliksty- 
Kochegov (BK) equation [223,224] that describes non-linear small-a; evolution in the 
saturation regime. This correspondence has been studied in detail in Refs [225,226], 
where BMS and BK were related via a stereographic projection. Because a 
generalization of the BK equation to hnite Nq exists [227,228], the correspondence 
between non-global logarithms and small-x physics was argued to hold at hnite-iVc and 
numerical solutions have been studied [229,230]. Very recently, this correspondence 
was indeed mathematically established [231]. We shall briehy discuss this topic in 
section 5.3. 


5.2. An alternative approach: soft-jet expansion 


A diherent approach to the question of resumming non-global logarithms was 
developed in Refs. [105,106] and applied to a phenomenological study of jet vetoes 
between hard jets in Refs. [107,232]. In that context, because color-correlations were 
of primary interest, the large-limit did not seem adequate. 

In this approach, the all-order calculation was then organized diherently, in terms 
of the number of soft (real or virtual) gluons in the region X, where the observable 
has no support: 


keX. 

nk=0 


(91) 


This procedure hxes the number of external partons (the soft gluons in X plus the 
original hard legs) and for each of these conhgurations an all-order evolution is derived. 
The kinematic structure of the latter is actually relatively simple, because after real- 
virtual cancellation, we are just left with virtual corrections in the X region, which 
do exponentiate. For = 0 we have no gluons in X, which corresponds to the 
resummation of the global component. Even if one is formally able to write all-order 
expressions for each of the in Eq. (91) in terms of abstract color operators, 

their actual evaluation requires hxing color bases and Ending representations for the 
color products. As already mentioned, this poses computational issues, especially 
because the dimensionality of the bases increases with the number of (real) gluons in 
X. Therefore, only the case n = 1 was considered in the phenomenological analyses 
of Refs. [107,232]: 
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(92) 


where the operator describes the change in color, spin and kinematics due to the 
emission of a soft gluon with momentum k. The first line correspond to a virtual 
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Figure 8. Pictorial representation of the color density matrix (j[U]. For each colored 
final state, an independent color rotation U°‘^{6) is applied between the amplitude 
An and its complex conjugate. Figure adapted from Ref. [231]. 


emission, which does not change the dimensionality of the color space. Therefore, the 
system evolves with the same resummed exponent as the global contribution 
where the superscript counts the number of gluons in X, zero in this case. However, 
in case of a real emission, as in the second line, we have now the evolution of a larger 
system, which is controlled by a different resummed exponent 

We hnish this discussion pointing out that an approach similar in spirit was 
recently developed in the context of jet substructure using techniques of SCET [233], 
where the equivalence between the soft-jet expansion (in the large-A^c limit) and 
an iterative solution of the BMS equation was established. Moreover, the soft-jet 
expansion was found to converge rapidly to the Monte Carlo solution of Ref. [71] . 

5.3. Recent developments 

In section 6 we will discuss the growing interest in jets physics and their substructure. 
This not only has lead to a rich phenomenological and experimental program in LHC 
physics, but it has also triggered a renaissance of all-order QCD calculations to better 
understand new ideas and techniques. The long-standing issue of resumming non- 
global logarithms has greatly benehted from this renewed interested. A remarkable 
study has been recently presented in Ref. [231]. In this approach, the color density 
matrix cr\U] is introduced, with the aim of describing soft radiation, as pictorially 
described in Fig. 8. An evolution equation (henceforth the Caron-Huot equation) is 
then derived for u\U], to all-loops, at hnite Nc- The related anomalous dimension 
K is explicitly computed to one and two loops. The one-loop approximation to the 
Caron-Huot equation coincides with the BMS equation, once the large-A'^c limit is 
taken and it confirms on a firmer ground the results of Refs. [229,230] at finite Nc- 
More importantly, the explicit calculation of the two-loop contribution to K paves 
the way for the resummation of non-global logarithms at higher-logarithmic accuracy, 
although computing solutions to the evolution equation remains a challenging task. 

6. Jet physics 

Jets, i.e. collimated sprays of particles, are key objects in particle physics. Jets really 
live at the boundary between experimental and theoretical particle physics and are 
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abundantly used by both communities. Indeed, the majority of physics analyses from 
the ATLAS and CMS collaborations uses jets as input. Consequently, jet dehnitions, 
commonly referred to as jet algorithms, have to make sense from an experimental 
viewpoint as well from a theoretical one. For instance, jet algorithms shonld be 
IRC safe, so that they yield hnite cross-sections when evaluated in perturbation 
theory [25] jj. 

Modern jet algorithms are based on the concept of sequential recombination. 
Pairwise distances between particles are evalnated in order to decide whether to 
recombine two particles. The metric used to evaluate these distances characterizes the 
jet algorithms. Nowadays, the most popular group of jet algorithm is the generalized 
kt family, for which the metric is defined by 

dij = min dw = Ptf, (93) 

where PuiPtj are the particles’ transverse momenta and Ai?^^- is their distance in the 
azimuth-rapidity plane. R is an external parameter, which plays the role of the jet 
radius. Different choices for the parameter p are possible. For instance, p = 0 
corresponds to the so-called Cambridge-Aachen (C/A) algorithm [235,236], with a 
pnrely geometrical distance. For p = 1 we have the /cj-algorithm [237,238] , which by 
clustering particles at low pt first, is likely to faithfnlly reconstruct a QCD branching 
history. Finally, with the choice p = —1 we obtain the dnii-kt algorithm [239], which 
clusters soft particles around a hard core, producing fairly round jets in the azimnth- 
rapidity plane. It is interesting to note that all algorithms of the generalized kt family 
act identically on a configuration with just two particles: they are recombined if 
ARij < R. 


6.1. Jet masses 

Many observables have been devised to study the internal properties of jets. Some of 
them were originally dehned to describe the properties of entire events, and then 
adapted to probe jets. Examples inclnde jet shapes [240], angnlarities [241,242], 
and energy-energy correlation fnnctions [243] of high-p^ jets. Perhaps the simplest 
example of such observables is the jet invariant mass. More precisely, we can define 
the dimensionless ratio 



The above ratio is small in the boosted regime mjet ^ Rpti which is of particnlar 
interested at the LHC, as discussed later in section 6.2. Thus, in order to obtain 
reliable predictions for this observable, we need to perform all-order calculations. 

The resnmmation of the p distribntion is clearly closely related to the resnmmation 
of thrust and of the hemisphere mass, which we have described in detail in the previous 
sections. Because large-angle radiation contribntes to single-logarithmic accuracy, in 

U For the interested Reader, we recommend the comprehensive review on jet physics of Ref. [234]. 
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order to achieve NLL, we have to take into account soft gluons emitted from each 
of the dipoles. However, in order to make our discussion simpler, we can work in 
the small-i? limit, neglecting contributions that vanish as powers of R. In this limit 
a simple picture emerges because large-angle radiation from dipoles other than the 
one involving the measured jet is suppressed [216]. Corrections to this picture can 
be systematically included [217] as a power series in the jet radius R. Moreover, 
logarithmic contributions at small-i? can also be resummed [244]. Thus, the NLL 
cumulative distribution for an isolated jet reads 


S(p) 


exp 




r(l-(2«,)-'ai(7i(Ap)) 






(95) 


where in case of a quark jets the functions Qi are dehned in Eq. (67) and the factor 
of 1/2 appears because thrust essentially corresponds to sum of the two hemisphere 
masses. 

We have already discussed the non-global factor and its resummation in 

section 5. Here we limit ourselves to note that non-global logarithms do not vanish 
in the small-/? limit. Moreover, the detailed form of the non-global contributions also 
depends on the clustering algorithm that dehnes the jet. In particular, in the presence 
of many soft emissions together with a hard parton, the anti-fc^ algorithm will always 
cluster all soft gluons to the hard parton, behaving as a rigid cone algorithm. In this 
case the algorithm provides a sharp boundary and non-global logarithms turn out to 
be the same as in the hemisphere case, up to corrections suppressed by powers of the 
jet radius. 

The function accounts for departure from the rigid-cone situation and 

therefore = i. The choice of different algorithms, such as C/A or kt can 

result into a reduction of the non-global contribution [245], however these algorithms 
are characterized by a non-trivial [216,246,247]. 

The jet mass distribution of the hardest jet in associated production of a Z boson 
with at least one jet, plotted as a function of C = R\/p is shown in Fig. 9. Jets are 
dehned using the anti-kt algorithm. 

The plot in Fig. 9, on the left, shows the NLL resummation in three different 
approximations. The result in the isolated jet approximation of Eq. (95) is plotted 
in blue. Corrections due to the hnite size of the jet radius are accounted for in the 
green curve. We note that all dipoles (three in this case) contribute now and the 
bulk of these corrections come from initial-state radiation. Finally, the resummation 
of non-global logarithms is accounted for in the red curve. 

The plot in Fig. 9, on the right, shows instead a comparison of the hxed-order 
0(as) result (black), the resummed one (green) and the matching between the two 
(red). At small values of (, the hxed-order curve exhibits the expected logarithmic 
divergence and in order to obtain a reliable prediction we have to resum these 
contributions to all orders. However, the resummed result is based on the eikonal 
approximation which is not accurate when mjet ~ Rpt- We can obtain a reliable 
description of the entire spectrum by matching the resummation to the hxed-order 
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Z+jet, R=0.6, pjj > 200 GeV 


0.15 

? = mj/pjj 



Figure 9. The jet mass distribution of the hardest jet in associated production of a 
Z boson with (at least) one jet, plotted as a function of C = R^/p- The plot on the 
left shows the NLL resummation in three different approximations: in the small-i? 
limit (blue curve, labelled as “Jet Functions”), with finite corrections (in green) 
and with non-global logarithms (in red). The plot on the right shows the fixed-order 
0{as) calculation (dotted black curve), the NLL resummation (green dotted curve) 
and the matching of the two (solid red curve). Both figures are taken from Ref. [217]. 


result, as shown in red. 

6.2. Grooming and tagging algorithms 

Because of its unprecedentedly high colliding energy, the LHC is reaching energies far 
above the electro-weak scale. Therefore, analyses and searching strategies developed 
for earlier colliders, in which electro-weak scale particles were produced with small 
velocity, have to be fundamentally reconsidered. In particular, in the context of jet- 
related studies, the large boost of electro-weak objects (not only Standard Model 
particles like electro-weak and Higgs bosons and top quarks, but also any new particle 
with a mass of the order of the electro-weak scale) causes their hadronic decays to 
become collimated inside a single big jet [248,249]. 

This is particularly important in the context of Higgs physics, because its 
dominant decay channel is into &-jets, which suffers from a huge QCD background. 
However, when the Higgs is produced with large transverse momentum, its decay 
products are likely to be reconstructed in one big jet. The presence of the Higgs 
can be then inferred by studying the substructure of this jet [250]. Consequently, jet 
substructure has emerged as an important tool for searches at the LHC and a vibrant 
held of theoretical and experimental research has developed in the past few years, 
producing a variety of studies and techniques [251-254]. 

Many “grooming” and “tagging” algorithms, e.g. the mass-drop tagger [250], 
trimming [255], pruning [256,257], have been developed, successfully tested and are 
currently used in experimental analyses. Broadly speaking, a grooming procedure takes 
a jet as an input and tries to clean it up by removing constituents which being at wide 
angle and relatively soft, are likely to come from contamination, such as the underlying 
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event or pile-up. A tagging procedure instead focuses on some kinematical variable 
that is able to distinguish signal from background, such as, for instance, the energy 
sharing between two prongs within the jet, and cuts on it. Many of the algorithms on 
the market usually perform both grooming and tagging and a clear distinction between 
the two is difficult. 

Regardless of their nature, substructure algorithms try to resolve jets on smaller 
angular and energy scales, thereby introducing new parameters. This challenges our 
ability of computing predictions in perturbative QCD. Indeed, most of the theoretical 
studies of substructure tools have been done using Monte Carlo parton showers. While 
these are powerful general purpose tools, their essentially numerical nature offers little 
insight into the results produced or their detailed and precise dependence on algorithm 
parameters. In the past few years however, we have reached a deeper understanding 
of, at least, the most basic (and most used) grooming and tagging techniques, both in 
the presence of background [258,259] and signal jets [260,261]. This understanding has 
been put at work and new substructure algorithms, which combine efficient signal / 
background discrimination together with robust theoretical understanding, have been 
devised. One of them is soft-drop [33], which we will use as an example for our 
discussion. 

The soft-drop procedure takes a C/A jet and implements the following steps: 

(i) Break the jet j into two subjets by undoing the last stage of C/A clustering. Label 
the resulting two subjets as ji and j 2 - 

(ii) If the subjets pass the soft drop condition > Zcut fhen deem j 

to be the final soft-drop jet. 

(hi) Otherwise, redefine j to be equal to subjet with larger pt and iterate the procedure, 
(iv) If j is a singleton and can no longer be declustered, then one can either remove j 
from consideration (“tagging mode”) or leave j as the final soft-drop jet (“grooming 
mode”). 

The difficulty posed by substructure algorithms in general, and soft drop in 
particular, is the presence of new parameters that slice the phase-space for soft gluon 
emission in a non-trivial way, resulting in potentially complicated all-order behavior 
of the observable at hand. This is exemplihed in Fig. 10, where we show the phase 
space for a single gluon emission from a hard quark/gluon in the (ln^,ln^) plane, 
where 0 < z < 1 is the energy fraction of the emitted gluon with respect to the hard 
parton initiating the jet and 0 < 6* < R is the angle of the emission, measured from 
the hard parton. This figure is analogue to the one we have used Fig. 5 to describe 
the allowed phase-space for event shapes, but now expressed in (ln^,ln^) variables 
rather than (In either set of coordinates the emission probability is flat in the 

soft-collinear limit. 

In the soft limit, the soft drop criterium reduces to 



^ In - < In-h /3 In —. 

2 : ^cut 0 


( 96 ) 
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log 



Figure 10. Phase space for emissions on the (lni,ln^) plane. For 13 > 0 , soft 
emissions are vetoed while much of the soft-collinear region is maintained. For 
13 = 0 , both soft and soft-collinear emissions are vetoed. For /3 < 0, all (two-prong) 
singularities are regulated by the soft drop procedure. Figure taken from Ref. [33]. 

Thus, vetoed emissions lie above a straight line of slope fi on the (ln^,ln^) plane, 
as shown in Fig. 10. For /3 > 0, collinear splittings always satisfy the soft drop 
condition, so a soft-drop jet still contains all of its collinear radiation. The amount of 
soft-collinear radiation that satisfies the soft drop condition depends on the relative 
scaling of the energy fraction 2 ; to the angle 6. As /3 —>■ 0, more of the soft-collinear 
radiation of the jet is removed, and in the fi = 0 limit, all soft-collinear radiation is 
removed jj. Therefore, we expect the coefficient of the double logarithms in observables 
like groomed jet mass, the origin of which is soft-collinear radiation, to be proportional 
to fi. In the strict /3 = 0 limit, collinear radiation is only maintained it z > Zcut- 
Because soft-collinear radiation is vetoed, the resulting jet mass distributions will only 
exhibit single logarithms, as emphasized in [258,259]. Moreover, non-global logarithms 
were found to be power-suppressed for /? > 0, and absent for /3 = 0. Finally, for /3 < 0, 
there are no logarithmic structures for observables like groomed jet mass at arbitrarily 
low values of the observable. For example, fi = —1 roughly corresponds to a cut on 
the relative transverse momentum of the two prongs under scrutiny. 

The above understanding can be formalized into actual calculations and the 
resummation of a variety of observables measured on soft-drop jets has been performed 
to NLL [33]. Examples of such observables include the jet mass, energy-correlation 
functions, the effective jet radius after grooming, the energy removed by the grooming 
procedure, as well as the energy sharing of hrst splitting that passes soft drop [34,262]. 
Two examples are shown in Fig. 11: on the left the groomed energy-correlation 
function, and on the right the groomed jet radius, that is the distance in azimuth 
and rapidity of the first pair of prongs that passes soft drop. The different curves are 
for different values of the soft-drop parameter fi, from no grooming (i.e. fi = 00 ), down 

ft Soft-drop with /3 = 0 corresponds to the modified Mass Drop Tagger [258,259]. 
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to very aggressive grooming (negative fi). 

One interesting feature that can be observed is the presence of a transition point 
in the distributions. Its position can be analytically predicted as a function of the 
algorithm’s parameters and jet’s kinematics. Moreover, the logarithmic structure of the 
distribution changes when we move across this transition point. In order to understand 
both of these features is actually enough to perform a very simple calculation. We 
consider the O (os) contribution to differential distribution for the energy-correlation 
[243] in both soft and collinear limits, i.e. z <^1 and 9 R\ 


1 da'"^ ‘^O-sCi 




(a) 



(97) 


where C* = Cf(Ca) in case of a quark (gluon) initiated jet. The S function hxes a 
measured value of C[°'\ in the soft-collinear limit, while the 0 function enforces the 
emission to pass soft drop, see Eq. (96). 

For /3 > 0, the evaluation of the two integrals is straightforward: 
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while for /3 < 0, there is an additional restriction which imposes a minimum allowed 
value for the observable Thus, we see that the result has indeed 

a transition point: for > Zcut the distribution is insensitive to the grooming 

procedure, while soft drop becomes active for < ^cut- Moreover, the coefficient 
of the double logarithm in this region is proportional to fi and consequently vanishes 
in the /3 —)■ 0 limit, as previously anticipated. It is also interesting to note that a 
direct O (ck^) shows that non-global logarithmic contributions to the spectrum 
are actually power-suppressed in the soft drop region < Zcut- 

The results discussed in this section, demonstrate that our understanding of 
the held of jet substructure has reach a level of maturity comparable to what has 
been achieved for more traditional observables, like the ones we have previously been 
discussing in this review. 


7. An unsettling end: breakdown of factorization 

The concept of factorization, i.e. the ability of separating physical effects that happens 
at different energy scale is the foundation of all the resummation program we have 
been discussing so far. Even more generally, we can say any QCD calculation, being it 
done at hxed-order or at the resummed level requires some notion of factorization. Of 
particular importance is the collinear factorization theorem [263-267] that allows us to 
separate the perturbative, i.e. calculable, part of a process from the non-perturbative 
one, which can be described in terms of parton distribution (or fragmentation) 
functions. These objects are universal, i.e. they do not depend on the particular 
process, and can be determined by htting data from previous experiments. 
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Analytic Analytic 



Figure 11. Two examples of NLL resummation for groomed observables. The 
grooming procedure of choice is soft-drop [33] . The plot on the left shows the energy- 
correlation function [243] after soft drop, for different values of f3. The plot on the 
right shows the groomed jet radius, i.e. the distance in azimuth and rapidity of the 
first pair of prongs that passes soft drop. Both plots are taken from Ref. [33]. 


Although explicitly proven only for a very few inclusive process, such as, 
for instance deep inelastic scattering of an electron off a proton [268], the Drell- 
Yan invariant mass and rapidity distribution [263,264,269], and heavy quarkonium 
production [270], collinear factorization is usually considered valid and is used 
ubiquitously in perturbative QCD calculations, regardless of their inclusive or exclusive 
nature. For instance, Monte Carlo parton showers are built on the idea of all-order 
factorization. 

The Authors of Ref. [271] embarked in an analysis of factorization properties of 
partonic scattering amplitudes and their related cross-sections. The analysis shows 
that strict collinear factorization of QCD amplitudes is violated beyond tree-level 
for initial-state (space-like) splittings. These factorization-breaking contributions 
originate from the exchange of Coulomb modes between the two incoming partons 
long before the hard interaction or between two outgoing partons, much later than 
the hard interaction [272]. As we have already seen, this results into an imaginary 
part, which cancels when evaluating the one-loop squared amplitude. However, as 
mentioned in section 4.2, this cancellation is not guaranteed at subsequent orders if 
the color-structure of the process is non-trivial. In particular, for the case of 2 —)• 2 
scattering of QCD massless partons, factorization-breaking contributions are found 
to cancel at two-loops, but they appear to give a non-vanishing contribution to the 
squared amplitude at three-loops and beyond. 

Let us then discuss what happens if we attempt to perform an all-order calculation 
of the (cumulative) distribution of an IRC safe observable v. Because we sum 
over all modes that live at a scale below v fully inclusively, the actual singularities 
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cancel between real and virtual corrections, ultimately because the colliding hadrons 
are colorless states. However, because we are measuring v, we have to veto real 
emissions above that scale. This, as discussed at length, induces potentially large 
logarithms, and prevents also potential factorization-breaking effects to cancel. The 
consequence is remarkable: the familiar picture of a hadron-hadron cross-section made 
of a hard scattering coefficient function and two universal PDFs, which are evolved 
independently up to the hard scale, breaks down (at high enough perturbative orders) 
if we insist on making the cross-section more exclusive by measuring an observable at 
a lower scale. 

The effect discussed above can have spectacular consequences on specihc 
observables. Let us consider, for instance, dijet production in hadron-hadron collisions. 
We identify the two leading jets and veto the presence of a third jet in the region (in 
azimuth and rapidity) between the leading ones. The cross-section for this observable is 
affected by large logarithms of the veto scale over the hard jet transverse momentum, 
which we wish to resum. The Authors of Refs. [105, 106] were set to tackle this 
project when they found an explicit example of factorization breaking. According to 
the standard analysis performed so far, cross-sections with central jet vetoes are non- 
global, because phase-space is restricted only between the hard jets. The cross-section 
is expected to be single-logarithmic, i.e. with all the logarithms originating 

from soft emission away from the jets. However, a super-leading contribution 
is found for n > 4, when the transverse momentum of the emission is chosen as 
the ordering variable [105,106]. The extra logarithms originate from the fact that 
emissions outside the region between the leading jets can go forward, i.e. collinear to 
the initial-state. However, these collinear contributions are prevented to cancel against 
the correspondent virtual contributions precisely because of the presence of Coulomb 
exchanges. 

Finally, we mention that the effect of factorization-breaking on global event shapes 
was discussed in Ref. [76], where it was noted that the logarithmic accuracy at which 
these effects kick in depends on the parameters hi (and possibly a) and crucially on 
the choice of the ordering variable. If transverse momentum was conhrmed to be the 
correct choice, then the effect would be very suppressed for 5; < 0 but could potentially 
be relevant, even at NLL, for hi > 0. 

8. Summary and Conclusions 

In this topical review we have discussed the basic concepts behind all-order calculations 
in QCD. We have focussed on calculations that make use of resummation in the so- 
called direct QCD framework. Similar results can also be obtained using the methods 
of Soft-Collinear Effective Theory. 

In the hrst part of this review, we have discussed basic properties of gauge- 
theory amplitudes in the soft and collinear limits, specihcally QED and QCD. We 
have started by reviewing in section 2 factorization properties at one loop, and then 
we have generalized this analysis to all orders in section 3. In the second part of this 
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review, we have instead focussed on phenomenological applications. In section 4, we 
have applied the all-order techniques previously presented to the event-shape variable 
thrust, discussing its NLL resummation in some detail. We have then presented the 
Caesar framework, which enables one to perform the resummation of a large class 
of event shapes and briefly mentioned its extension to NNLL. Section 5 was instead 
devoted to a discussion of non-global observables and their resummation. In section 6, 
we have illustrated how logarithmic resummation is an invaluable tool to describe 
the physics of hadronic jets and, in particular, their internal structure, which is a 
topic of the highest importance in the context of new physics searches at high-energy 
colliders. Finally, we have briefly discussed important issues about the limitations of 
the framework we have used to factorize amplitudes and cross-sections. 

With this review, we hope we have managed to stimulate the curiosity of a Reader 
who begins to study this field, explaining them the concepts and the methods on 
which all-order calculations are based, with a constant eye on their phenomenological 
implications, as well as recent developments in the held. We realize that a full account 
of technical details would require a text-book, rather than a review. Therefore, 
we have put an effort in compiling an exhaustive bibliography and we encourage 
the Reader to read the original papers in order to gain a deeper, more technical 
understanding. Because we feel that a concise summary of this kind of calculations 
and their applications was lacking, we also reckon this effort may prove useful also for 
a more expert Reader. 
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